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ABSTRACT 


The  complete  and  rigorous  update  procedures  applied  to  the  1950.0 
epoch  SAO  Catalog  are  described.  The  corrections  included  are  the  terms 
of  elliptic  aberration,  precession,  proper  motion,  and  foreshortening 
(e.g.,  parallax  and  radial  velocity  terms).  The  uses  of  the  updated  catalog 
are  schematically  described  for  the  local  reduction  of  artificial  satellite 
positions.  Finally,  an  annotated  FORTRAN  program  is  included. 


iii 


I.  INTRODUCTION 

This  report  describes  the  sequence  of  computations  performed  to  update 

3 

the  Smithsonian  Astrophysical  Observatory  Catalog  (SAO  ; hereinafter  the 

SAOC)  from  epoch  1950.0  to  the  present.  We  also  schematically  describe  the 

use  of  a current  version  of  the  SAOC  in  reducing  artificial  satellite 

4 

positions  and  photometry.  The  reader  should  be  familiar  with  Taff  . Also, 

Taff  (in  press^  discusses  photometric  reductions.  Finally,  the  reader  should 

4 

be  aware  that  the  errors  discussed  in  § VIII  of  Taff  are  reduced  by  a factor 
of  two  when  proper  treatment  of  the  e-terms  is  included. 

Sections  III  and  IV  detail  the  update  calculations  performed  on  the 
SAOC  and  Section  V describes  the  astrometric  information  contained  in  the 
1975.0  master  and  subsequent  yearly  tapes. 
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II.  DESCRIPTION  OF  THE  SAOC 


The  SAO,  in  response  to  the  needs  of  the  space  age,  constructed  a non- 
fundamental, astrometric,  compilation  catalog  with  the  following  goals: 

(1)  the  star  density  would  be  at  least  four  stars  per  square  degree,  (2)  the 
1970.0  positions  would  have  a root  mean  square  error  of  £ 1",  and  (3)  the 

catalog  would  be  on  the  FK4  system.  The  physical  existence  of  the  SAOC  and 
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reference  to  Haramundanis  and  Scott  and  Smith  indicate  that  these 
goals  have  been  met. 

The  catalog  contains  entries  for  258,997  stars.  For  astrometric  purposes 
the  most  important  aspects  of  the  SAOC  are  its  success  in  meeting  the  above 
goals  and  its  availability  in  machine  readable  form.  For  our  purposes  the  most 
important  datum  for  each  star  consists  of  position  (ot  * right  ascension, 

6 = declination)  and  proper  motions  (y  - proper  motion  in  right  ascension, 
y1  = proper  motion  in  declination).  This  information  is  in  a barycentric 
solar  system  coordinate  system  for  the  equator  and  equinox  of  1950.0.  The 
reader  is  referred  to  the  Introduction  of  the  SAOC.  However,  he  should  be 
aware  that  the  column  is  totally  incorrect  and  although  the  wrong  equations 

appear  on  p.  xiii  for  cr^  and  , , the  correct  formulas  were  used  in  their 
computation . 
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III.  ELLIPTIC  ABERRATION 


The  position  (e.g.,  a,  6)  of  a star,  referred  to  a mean  equator  and 
equinox  of  date,  is  referred  to  as  a mean  position.  It  differs  from  the 
true  position  (at  that  time  in  that  coordinate  system)  by  the  effects  of 
nutation.  However,  because  the  eccentricity  of  the  earth's  orbit  is  small 
and  hand  computation  laborious,  it  has  been  traditional  to  leave  in  those 
observational  effects  due  solely  to  the  ellipticity  of  the  earth's  orbit  in 
the  just  defined  mean  place.  Catalogs  contain  mean  places  of  stars.  Let 
[aQ(t ' ) , 6Q(t')]  be  the  just  defined  mean  place  at  epoch  t'  as  found  in  a 
catalog.  Then  the  real  (i.e.,  intrinsic)  position  [a(t'),  6 ( t ' ) ] is  obtained 
by  setting 


a(t ' ) * a (t ' ) - Aa(t ' ) , 
o 

6(t ' ) = 6 (t?)  - A6(t ' ) . 

o 


(la) 


(lb) 


The  e-term  corrections  are  given  by 


Aa(t ' ) * cAC  + dAD , 


(2a) 


A6(t ' ) - c'AC  + d'AD, 


(2b) 


where  the  Besselian  star  constants  c,  d,  c',  d'  are  given  by 


c = cosa  (t')sec6  (t'),  d = sina  (t')sec6  (t'), 
o o o o 

c'  = tanc(t')cos6  (t')  - sina  (t')sind  (t'), 
o o o 


(3a) 


d'  = cosa  (t')sin6  (t'). 


(3b) 


o 


o 


3 


where  e is  the  true  obliquity  of  the  ecliptic  at  t * t*.  The  quantities 
AC  and  AD  are  computed  from  e,  the  constant  of  aberation  K(=  2074958) , the 
eccentricity  of  the  earth* s orbit  e,  and  the  longitude  of  the  earth *s 
perihelion  0).  Thus, 

AC  = e ( t * )kcosu)(  t * )cose  ( t * ) 

AD  = e(  t ' )Ksino)(  t * ) , 

where 

e = 0 . 01675104  - A. 18  x 10_5T  - 1.26  x 10_7T2, 
e = 23°27 ' 8'. '26  - 467845T  - 070059T2  + 0V00181T3, 
to  » 101°13' 1570  + 6189703T  + 1763T2  + 07012T3, 

and  T is  the  number  of  tropical  centuries  elapsed  since  1900.0. 

From  Eqs . (5)  it  is  clear  that  e,  e,  and  U)  vary  very  slowly  and  that 

eK  is  a small  number.  Hence,  these  corrections  are  generally  minute.  In 

fact,  it  can  be  shown  that  if  t’  35  1950.0,  tM  = 1975.0  then  the  absolute 

value  of  the  maximum  errors  in  a(t")  and  5(t")  due  to  neglect  of  the  e-terms 
-3  2 -3 

are  2711  x 10  sec  6(t’)  and  0797  x 10  sec6(t*).  The  errors  are  linear  in 

t"  - t*. 

However,  as  is  the  case  with  most  approximations  used  in  spherical  astronomy, 
the  neglect  of  the  e-terms  produces  an  effect  in  right  ascension  that  systematically 
increases  with  declination.  Therefore,  since  these  right  ascensions  are 
necessarily  of  poorer  accuracy  than  those  at  lower  declinations,  the  first 


(4a) 

(4b) 


(5a) 

(5b) 

(5c) 


4 


step  in  updating  the  SAOC  is  to  apply  Eqs.  (1).  The  second  step  is  discussed 
below  (§  IV).  The  third  step  is  to  add  back  in  Aa(t"),  A6(t"). 
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IV.  PRECESSION  AND  PROPER  MOTION 


Using  a(t')  and  6(t')  [along  with  y(t')  and  y'Ct1)]  as  input,  the  power 

4 

series  method  of  Taff  was  then  applied  from  tf  = 1950.0  to  t"  = 1975.0. 

The  complete  details  are  in  § III  of  that  report.  The  Appendix  contains  an 
annotated  listing  of  a thoroughly  tested  FORTRAN  program  to  perform  all  of 
the  computations  discussed  herein. 

4 

While  only  minor  typographical  errors  have  been  discovered  in  Taff  , the 
algebraic  reduction  which  yielded  Eqs.  (14a,  b)  is  sufficiently  complicated 
that  the  unsimplified  versions  are  used  in  the  program. 
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V. 


THE  LINCOLN  1975.0  YEARLY  VERSIONS  OF  THE  SAOC 


After  a complete  reduction  had  been  performed  to  1975.0  a new  SAOC 
tape  was  constructed  at  Lincoln.  It  contains  for  epoch  1975,  a,  5,  y,  y’, 
da/dt,  0.5d^a/dt^,  d6/dt,  0.5d^6/dt^,  dy/dt,  and  dy'/dt.  In  addition, 
depending  on  6 ( t f ) (relative  to  ± 80°),  either  1950.0  0.5d^a/dt^,  0.5d^6/dt^, 
dy/dt,  and  dyf/dt  or  the  same  quantities  at  1970.0.  The  tape  has  been  resorted 
into  10°  wide  declination  bands  with  right  ascension  increasing  in  each  band. 

For  each  year  the  1975.0  tape  is  updated  (on  or  about  1 July)  to  provide 
positions  and  proper  motions  whose  epoch  is  that  of  nearest  Besselian  solar 
year.  This  update  is  performed  using  Eqs.  (18,  19)  or  (20,  21)  of  Taff^. 
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VI.  USES  OF  THE  SAOC 


The  primary  use  of  the  SAOC  is  in  the  precise  reduction  of  artificial 
satellite  positions  and  photometry.  For  this  purpose  an  interactive  graphics 
display  program  exists.  The  operator  (or  computer)  enters  the  right  ascension 
and  declination  of  interest,  the  size  of  the  surrounding  area  to  be  displayed 
(up  to  5°),  and  either  the  star  density  per  square  degree  (with  the  brightest 
stars  used  first)  or  a magnitude  limit.  The  program  then  generates  a display 
of  the  appropriate  star  field  including  differing  intensities  for  the  stars. 

The  position  of  the  telescope  and  the  satellite  in  question  are  also  indicated. 
As  the  real  time  pointing  system  provides  information  on  telescope  motion, 
it  is  a simple  task  to  continually  update  the  display,  look  backward  in  time, 
or  look  forward  in  time.  A local  calibration  procedure  will  be  used  to 
simultaneously  reduce  the  position  and  photometry  of  the  satellite. 
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APPENDIX:  LISTING  OF  FORTRAN  PROGRAM 


Following  this  text  is  a complete  listing  of  the  two  subroutines  used  to 
update  the  SAOC.  In  particular  the  epochs  1950.0,  1975.0,  and  1977.0  appear 
explicitly  in  the  program  (e.g.,  within  DO  loops  //l,  it 2;  above  statement 
it 4;  the  IF  statement  below  statement  it 5;  at  statement  it 6;  at  statement  //14;  and 
below  statement  // 15 ) . Quantities  such  as  R(I)  , RA,  ALP  refer  to  right  ascension; 
D (I) , DEC,  DEL  refer  to  declination;  MU,  MUA  refer  to  proper  motion  in  right 
ascension;  MUP , MUD  refer  to  proper  motion  in  declination;  DA,  EA  to  da/dt,  etc. 
Trigonometric  functions  of  an  angle  are  usually  prefixed,  thus  CO  is  the  cosine 
of  omega,  TD  the  target  of  declination,  etc.  With  close  reference  to  Eqs . (7-21) 

4 

of  Taff  the  reader  should  have  no  difficulty  following  the  program. 
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SUBROUTINE  uPUATE  <Rd ) tOCl ) tHUd ) tHUPd )«RC6)*D(6)«MU(6) • MUP<  6) ) 

C RIGOROUS  PRECESSION  $ PROPER  MOTION  UPDATE/POWER  SERIES  METHOD 
C INPLIf  ARE  VAKlARLtS  AT  INDEX  It  OUTPUT  ARt  VARIABLES  AT  INDEX  6 
C PRECISION  NEEDED  IS  10. **(-9)  FOR  ALL  OPERATIONS 
IMPLICIT  DOUBLE  PRECISION  (A-HtM-2) 

DOUBLE  PRECISION  KAP 
INTEGER  HRtUEbthINTtMlNA 

DIMENSION  R(fe>  tUA(b)  *U2A(6)  «U(6) tDD(6> ♦ 020(6) tMU(fe)  tUMU(fa)  tMUP<6) . 

* DMUP  ( 6 ) irt(b)  * MP ( 6 ) tN(6)  *NP(6)  * TE ( 2 ) t DELC ( 2 ) t UELD  ( 2 ) tPHOTdl  ) 

COMMON  /DERIV/UA7S.D2A50  t02A /5 * DU75 * D2U50 t D2D75 t DMU50 t DMU75 t 

* DruPbO tLMUP7b 

CDhM0N/UAT/P12tPI .TPI * T • F T * T WT 
C SET  UP  CONSTANTS  , 

DATA  KAptCU;Jt  iPtIM/2.049btiUltl.02269D-4f  ♦♦  •**-  •/ 

Pl  = 4. l41b<*2s5bbU0 
P 12  = 5 # U • 1 ♦P I 
1 Pl=Pl+Pl 

LKA  = Pl/4. 421)4 
LUEC=PI/b.460b 
PI  12  = P I / 1 • 2l  1 1 
PI180=Pl/l .602 
US=1 .DU/b.Dl 
I w r=2.ci/a.uo 
H=4.U0/3.D0 
UMrtX=« ,0D1*PI ItiO 

C I F PARALLAX  ANU  RADIAL  VELOCITY  UNKNOWN  SET  THEM  EUUAL  TO  2ER0 
PAK=U .00 
VK=0 .DO 

C SET  DP  PRECESSION  CUNSTANTS 
DO  1 I = 1 i b 
T=b.U-l+b.0-2*(  1-1  ) 

M(I)  =C  RAM3. 0/244/ U2M*(1. 8b  3UD-1  + T*  8. 0-6)  ) 
np(l)=CRA*(i.tibA0U-l+T*l.feU-5) 

N ( 1 )=CUEC*  <2. 0 0 468 5U3-T *(8.5330-14 T* 3. 70-4) ) 

1 ivp(l  )=-CUEC*(8.b34U-l4j  W.4U-4) 

MPP=CRA*1 .bO-b 
NPP=-CUEC*7.4U-4 

L SET  DP  tLLlPTIC  ABERRATION  CUNSTANTS 
DU  2 1=1*2 
r=5.U-l  + 2.50-l*d-l) 

UML=P1 18U*  ( 1 .01U2  + 0S*  ( 1 .3U140S*  ( 1.501  + T*(6.189Q3D3+T*(  1 . 6300  + T * 1 . 2 

♦ U-2  ) ) ))  ) 

CO  = LiCUS  (OML  ) 

SU=l'SlMUnE) 

EHS  = PU60*  (2.4U14US*  (2.7Ul40S*(8.26DO-T*(H.684  5Ul  + T*(5.9U-3-T*1.61 

♦ U-4) ) ) ) ) 

LE=OCOS(EPS) 

TE  (I  ) =0  T AN ( EPS ) 

LCC  = l .6  751040-2-T* (4 . 1 80- b+ T * 1 . 260 - 7 ) 

DELL ( 1 ) =LUEC*KAP*ELC*CO*Lt 

2 UELDt 1 ) =CUEC*KAP*ECC*SO 

C INPUT  IS  HA  ANU  DEC  IN  RADIANS?  MUtMUP  IN  RAO/CLNT  ? PAR  IN  SECONDS 

C OF  ARC*  VR  IN  K'VSEC 
L ULTERP1INE  OELLI^'AdON  RANGE 
L = b 

IMDAbS(Ud)  ) .GT  .UMAX)  L=1 
C ADJUST  KA  AND  DEC  FUR  E - TERMS 
5R  = LSIMR<  1 ) ) 

LR  = DCGS(Rd  ) ) 
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SU=LS1N  ( L>  ( X ) ) 

LU=CCOS  (1)11)) 

KI1I:KI1 ) -(CR*UELt ( 1 ) +SR*OELU(  1 ) )/CD 

uai:U|ll-(  ( Tt  ( 1 ) *CD-Sl<*SOI  *UEl_C  < 1 ) +CR*SD*DELD  ( 1 > ) 

T=5.U-2*L 

C INITIALIZE  COUNTER 
* = 1 

C CALCULATE.  UHA/UT  »U2Ra/DT**2,U3RA/0T**3,0DEC/UT,D2UEC/UT**2 ,D30EC/UT**3 
L » UMU/U 1 «u2ivU/UT**2.  UHUP/UT  • AND  U2MUP/DT**2  AT  1950.0 
A 5R  = DS1MR(K)  ) 

LK-llCOS  ( K (K  ) ) 

!U=DTAN(0(K) ) 

SEC2=1 .UO+TU*TU 
S2U=USIN(2.00*U(K)  ) 

C2U=UCoS(2.00*0(K)  ) 

PK=M(K)+N(K)*SK*TU 
PU  = N ( K ) *CK 
UA  (K  ) =PR+MU(R  ) 

UU ( K ) =PU+MUP ( K ) 

UPH=hF(K>+NP(K)*SK*T'J*N(K)*DA(K)*CR*TU+N(K)*00(K)*SR*SEC2 
UPU=NP  C K ) *Lr-N( K ) *0a  (K  ) *SR 

U2A(K)=b.0-l*(MP(K) + NP (K  > *SR*TO+N(K )*( DA (K )+MU( K ) ) *CR*  IU*N ( K ) * ( 00 ( 
*K ) +MUP (K ) )*SR*SEC2+2.UO*MU(K)*MliP(K)*TD-2.UO*MU(K)*CON*PAR*VR) 
U2UO\)=b.U-l*(NP(*)*CR-N(K)*(OA(K)+MU(K) )*SR-5.U-1*MU(K)*MU(K)*S20 
*-2  .UO*l-UPl  K ) *LUN*PAR*VK  ) 

UiILM  k ) =N(K )*MU(K)*CR*TO  + N(K)*MUP(K)* SR *SEC 2 + 2.00 *Mu(K)*MUP(K)*TD-2 
*.U0*NU(K  ) *CON*PAR*VR 

Ur, UP <K)=-N(*>*nU(K)*SR-5.U-l*MU(K>*MU(K) *S2U-2 . UO *MUP ( K ) *CON*PAR*V 

* K 

U3 A = MPP  + NPP * SR *TU  + NP(K ) *UA (K ) *CR*TU+NP ( K ) *00 ( K )*SR*StC2+Np (K ) * (PK+ 
*2.U0*Mu(k)  )+CK*lU*N(K)*(UPR+2.D0*DMU(K)  )*CR*TU”N(K)*(pR+2.00*MU<k) 
-o  )*UA(K)*SR*TL*N(K)*(PK+2.O0*nU(K))*OD(K)  *CR*SEC2+2  . UO  *OHU  ( K ) *hUP(  K 
ai)*TO  + 2.UO*Ml'(K)*DMUP(K)*TU  + 2»UO*MU(K)*MUP(K)*UO(K)  *SEC2+NP  ( K ) * < PD  + 
ri2.U0*MUP(K)  )*SH*SLC2*N(K)*(OPU+2.UO*OMUP(K>  ) *SR*SEC2  + M K ) * ( PO+2  . UO 
1 *MUP( K ) ) * U A ( K ) *CK* SEC 2+2.00 *N(k )*(PD+ 2.00 *MUP(K ) ) *SH*T0*SEC2*UU(K) 
u3U=NPP*CR-f;P(K)*UA(K)*SK-NP(K)*(PR+2.D0*UMU(K>  ) *SR-N ( K ) * ( UPK+2 . 00 

♦ *OriU  (M)*SR-MKl*l  PR  + 2 . UO  *MU  ( K ) ) *0A(  K ) *CR-MU (K  ) *OMU(K  ) *S2U-MU(K ) *M 
,?U(K)*UU(k)*c2U 

U2MU  = i*UI(  K ) *NP(K)  *LR*  TU  + MUP  ( K ) *NP(K ) *SR*SEC2+N(K )*(MU(K )*UU(K  )+HUP( 
*K  ) * 0 A ( K ) ) *CR*5EC2+n (K ) *(2.U0*MUP(K ) *00 (K ) *SEC2-MU(K ) *DA(K ) )*SH* TU- 
VA 2.UU*SU*SU*.v;U(K)**3+N(K)*N(K)*rtU(K)*(UR*CR*T0*TU-SR*SR*SEC2)*N(K)* 
o N ( K ) *M0P ( K ) *SH*CR*TU*SEC2  + H ,UO*N(  K ) *MU(K ) *MUP ( K ) *CR*TO*TU+N ( K ) * ( 2 . 
fjuO*MUP  ( K ) *Mi)F  ( K ) *SEC2-3.U0*MU(  K ) *MU(  K ) ) *SR*TU+4 . 00  *HU  ( K ) *MUP(  K ) *MU 
iP(R)*  fO*lU  + 2.UO*MU(K)*nuP(K)*00(K)*SEC2 
U2riuP  = -Mu  (K  ) *NP(K  ) *SR-N(K)  *MU(  K ) *DA  ( K ) *CR“N ( K ) *N  ( K ) *ML'P  ( K ) *SK*SR*S 
*EC2-Mk)*Mk)*MU(K)*SK*CK*TU-2.U0*N(K)*MU(M*HU(K)*CR*5>O*SO-‘*.U0*N 
ii  ( R ) *f*l)  (K  ) ♦ MUP  (K)*SR*TO-MU(K)*MU(K)*(  OU  (K  )+2«U0*(MUP(K  ) -PU ) *SU*SO) 

L OBTAIN  LSI  IMA ( LS  FUR  K A . UEC . MU , MUP  AT  1950. 0 + T 
b KA=P(K)+I*(UA(K)+I*(02A(K)+T*03A/6.DO)  ) 

ULC=U(k)+T*(DU(K) +T* (0201 K ) +T*03U/fe.OO ) ) 

HUA=nU(K )+T*( UMU(K )+5.0-l*T*U2MU> 

MUU=MUp ( K ) ♦ r * ( UMUP ( K ) *5 . 0-1* T * J2MUP ) 

SK  = DS1MKA) 

LR=ULOS(KA) 

( U=0T AN ( ULC ) 

SLL2=1 .UO  + TO*  TU 
S2U=OSIN(2.UO*ULC) 

C IF  T =0.23  UUNL 

lF(L.Ej.O)  GO  TO  19 

L T = O.Ob:  US  1 No  1955.0.1960.0,  ETC.  ESTIMATES  FOR  HA«UEC.MU,  ANU  MUP  RE- 
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C CALCULATE.  02RA/0T**2,U20EC/0T**2.0mU/0T.  ANO  OMUP/DT 
PK=M(K*1 )*N(.K*l)*SR»TO 
PU  = N ( K-fl  ) *CR 

L2A=5.0“l*(MP(K+l)+NP(K  + l)*SK»Tn+  N ( K + l ) * ( PR-t-2 . l)0*MUA ) *CK*TO  + N ( K + 1 
*)»(PU-*-2.U0*rtUU)»SK*SEC2  + 2.00*MUA*MU0*T0-2.U0*MUA*C0N*PAR*VR) 

L2U  = 5.0~1*  (NP(  K-t-1  ) »CR-  N(K*1  )*(  PR-t-2. 00*MUA)*SR -5.0-1  *MUA*MUA*S20-2 
*.DU*MUC*C0IJ*PAK*VK) 

UMUA  = N(  K-t-1 ) ♦r1UA*CH*TU-*-N(K  + l )*MuO*SR*SEC2+2.D0*MUa*MUO*TU-2.O0*MuU* 
*LON*PAR*VK 

UriUU  = -rj  ( K*  1 ) *MUA*SR-5.U-1*MUA*MUA*S20-2.U0*MUU*C0N*PAR*VH 
C GET  1955 . U , 19b0 . 0 « LTC.  HA. DEC. MU.  AND  MUP 

R<K  + 1 )=R(KH-T*(UA(K)-t-T*(U2A(KH-TWT*T*<E2A-02A<K)  ) ) ) 

U(K  + 1 )=U(K)+T*(UD(K)-t-T*<U20<K)-t-TWT*T*<E2U-02U(K)  ) ) ) 
rtU<K-H)=rlU<K)*T*U)MU<K)-l-l.Ul*T*<OMUA-UMU<K)  ) ) 

MUP  ( K-f 1 ) = MUP  ( K ) + T*  ( UMUP  ( K ) +1 . ul *T*  ( OMUO-OMUP  ( K ) ) ) 

C INCREMENT  COUNTER 
K = K + 1 
C UONE? 

b IF(K.Eu.b)  (-.0  IU  7 
C MOT  UUNL « REPEAT  UPUATE  CYCLE 
GU  TO  4 

C PRECESslUN  ANU  PROPER  MOTION  UPDATE  IS  NOW  UONE.  UNUO  E-TERMS 
7 SK=|)SlN<H<6) ) 

CK=DCuS ( R ( b ) ) 

SU=USlN(U(b) ) 

CU=UCOS ( U ( b ) ) 

R(b)=R  <b)F<CR*UELC  (2  ) +SR*OELU  ( ^ ) )/C0 

U(b  )=U(b  ) ♦<  ( TE(  2)  *Cu-SK*SU)*ULLC  (2)  -t-CR*S0*UELU(2  ) ) 

C CALCULATE  02RA/l)T**2,020EC/0T**2.0hU/0T.  ANO  DMUP/OT  AT  1975.0  FOR  USE 
C IN  YEARLY  IJPUAIL  PROGRAM  ( TO  YR ) 

SR  = OSI.\(K(K)  I 
LR=DCOS<K<K) > 

IU  = UTAIj(u(K)  ) 

StC2=l .UU-t-TU* IU 
S2U=USlN(2.uO*u(K) ) 

PR=M  ( K ) +1.  (K  I *SK*TU 
PU=N ( K ) *LR 
UA(k)=PR+Mu(K) 

UU(K)=PU+MuP(K) 

U2A<H)=5.U-l*<MP<H)+NP<K>*SR*T0-t-N(K)*<0A<K)-t-MU<K)  ) *CR*  TO-*-N  ( K ) * < OU  ( 
*K ) *MUP(K ) ) *SR*SEC2+2.00*MU(K ) *MUP(K) *TO-2.U0*MU(K )*CON*PAK*VK) 

U2U( h lsb.U-1* INP<K  »*CR-N»K»  *(OA(K» *MU<K ) ) *SR-5.U-1*MU ( K ) *MU(K ) *S2D 
*-2.U0*|4UP  (K  ) *CUN*PaR*VH  ) 

UPiU(K)=N(K)*MU(K)*Ck*T0*n(k)*MUP(K)*SK*SEC2+2.00*MU(K)*MUP(k)*T0-2 
* .UO*MU(K )»CuN*PAR»vR 

UMUP  <K )=-Nt  k )*MU(K ) *SR-5.U-1*MU(K )*MU«K) *S2D-2.U0*MUP (K )*CON*PAR*V 

*R 

C CHECK  FUR  OVER  THE  POLES  OR  ARUUNU  THE  VERNAL  EQUINOX 
lMUAbSlU(b)  ) .LT.P12)  GO  TO  11 
R(b)=R(6)*Pl 

IFCDCb) » tt.a.9 
a u(b)=-u(b»-pi 

GO  TO  10 
9 U(b)=PJ-U(b) 

10  IF t R (b  ) .GE . TPI ) R(b)=H<6)-TPl 
IFlK(b).LT.O.UO)  R ( fa  ) =H  ( b ) -F  TP  I 

L (JUT  PUI 

11  R A— R ( b ) /PI 12 
HR  = RA 

RA=b.Ol*(RA-HR) 
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J*1 1 NT  = K A 

KA=b.Ul*(HA-MlNf ) 

ULC  = UABS(U(6)  )/PH80 
UtbrOtC 

UtCrb.Ol* (UtC-ULG) 

M1NA=ULC 

LLC  = 6.iU*(UEC-MImA) 
nUArHU ( b ) /CRA 
MUU  = r.up  <fe  ) /CUEC 
1S=1P 

If-  (Li(b)  .LT.O.UO)  IS  = Irt 

WHITE (3 • 12)  HK •MINT ,KA,MUA, I S « OEG ♦ MI NA ♦ OEC « MUU 
12  h UKpi  A T UXi2l3tFtt.4i2XtF8.3i5XiAliI2iI3iF7.2«2XtF7,2) 
bO  Tu  lb 

C uSlrjb  l<*7b.U  ESTIMAIES  FOR  RA«OEC»Mu«  AMO  MUP  RECALCULATE  U2KA/01»*2« 

C f 2ULC/UT**2tU^U/0T • AND  DMUP/OT 
14  K = D 

PK  = PMK  + 1)+N(K  + 1)*SH*TU 
PU=Mh  + l )*CR 

L2A  = b.0-l*  ( .op ( K + l ) + ,\jP  ( k + 1 ) *SK*  TO  + MK  + 1 ) * ( PR  + 2 *U0*hUA  ) *CR* TO  + N ( K ♦ 1 
*)*(PU  + 2«U0*n:uU)  *SH*SEC2+2.UU*MUA*MUO*TU-2*OU*hUA*CON*PAK*VR  ) 
t2D  = b.»j-l*  <uP  (K  + l ) *CK-  M K + l > *<  PR*2.DU*MUA  ) *SR-D . 0- 1 *hU A*MUA*S20"2 
*,Oa*Muo*LUN»PAR*VK ) 

UttUA=M  K + l ) ♦MUA*CN*TU+N(K+1 ) *MUU*SR *StC2  + 2 . DO*MUA*MUD*  TU-2  • UU*KUU* 
•LON*PAR*VH 

UhuU  = -M  K + l ) *MUA#SR-b.0-l*riUA*f1uA*S2D-2.DU*MUU*C0N*PAR*VK 
C GET  lWb.O  RA«UEC«rtU«  ANU  rtUP 

h (6  >=k  < 1 ) +1  * IUA  ii  l-fl  * CU2A(  1 ) ♦FT*T*  (E2A-D2A  ( 1 ) ) ) ) 
U<b)=U(l)+T*<U0m+T*<02U(l)+FT*T*<E20-02UU)  ) ) f 
r.U(  b ) =hU(  1 ) + T*  < UfnU  ( 1 ) +2.U0*T*  < OMUA-OMU(  1 ) ) ) 
hUP( b ) =huP( 1 ) +1  * <UrtUP< 1 ) +2.U0*T* (OMUU-DMUP( 1 ) ) ) 

L St  T CUUMtrt  TU  END 
r\  = b 

bU  TO  t, 

ID  CONTINUE 
L TOYP  St L T 1 ON 

L MtKL  UPDATES  1 ^ 7b • U TU  1^77.0 
f =2.U-2 
U A 7 b - U A < b ) 
u2A7b  = l)2A  ( b ) 

DU7b=UU(fe ) 

U2U7b=J2U ( b ) 
unu7b  = 'JhU(b  ) 

U*UP7b=UhUP ( b ) 

L = b-L 

U2AbU=|)2rt  ( L ) 

U2UbU=02U (L ) 

U(-IUbU  = L;MU  (L  ) 

UhuPbUrO^OPd  ) 

CALL  T0TK(R(b)«U(b)«r1U(6)  «MUP<6)  « R A • OEC  ♦ MU  A ♦ MUD  « L ) 

C OUTPUT 

ALP=RA/P1 12 
MK=ALP 

ALP  = b.ljl1M  ALP“HK) 
niNT=ALP 

ALP  = b.Ul*t  ALP-Mit\lT  ) 

UtL=UAbS(ULC )/PllbO 
ULb=Ut L 

DEL=b.Dl*<ULL-UtG> 

rtlNA=OEL 
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L)EL  = 6 • 01  * ( UEL-M1NA  ) 
nUArnUA/CRA 

riUUrriUD/CULC 

1S  = IP 

IF(UEC.LT.O.UU)  IS=IH 

WRITE  ( 4 « ) HR  • MINT  « ALP1HUA1  IS ♦ DEG ♦ M I NA » DEL ♦ hUD 
16  CONTINUE 

call  exit 
Enl) 
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SUbPOUTlNL  TGYK ( ALP  tULL  tMUA  tMUOtHA  «DEC  «MU«MUP«L ) 
ir* PL1C1T  UOUbLL  PRECISION  (A-H»IW) 

LUiiF.ON  /UEK  IV /UA7btU2 AbO  tU2  A 75,007b  , 02050  • U2U7b  » OrlUbO  t 0MU7b  ♦ 
♦ unuPbu  »unup7b 

CUmhON/UAl/P12tPI « TPI * T «F I • TWT 
C INPUT  lb  ALP*ULL  IN  KAO;  MUAtrtUU  IN  RAD/CLNTI  T IN  CENT 
L CUMUl  IS  HAfCECthU»nUH  IN  SAME  UNITS 
L L IS  Mifeh  UtLLlNATlON  FLAG 
IF  (L.L(j.b)  GO  TU  1 

KA  = ALP  + T*  <UA7b+T* (U2A7b  + FT*T*( D2A7b-02AbO  ) ) ) 

ULC=DLl+T*  <UU7b  + TM02U75  + FT*T*(02U7b-U20bO ) ) ) 

MU  = hUA  + T*<  UNU/b«*2.*T*  (UMU7b-UMUbO  ) ) 
moP  = nUu+  I * ( Ol*Up75+2  • *T  * ( UfiUP  ?b-0«UP50  > ) 
b(J  TU  2 

1 HA=ALH+l*(UA7b+l*(O2A7b+lwT*J*(D2A75-D2Ab0) ) ) 

ULL=ULl+T  * (UU/b  + T* <Q2U7b+TWT*T* ( U207b-U20b0 ) ) ) 

PU=RUA+1  * ( DRU7b*  1U  . *T*  (OriU7b-OMUbO  ) ) 
nuPrRU^+7  * <MHUP7b+10. *T*  ( UhUP7 b-OMUPbO  ) ) 

C Ofc'C*  FCK  OVER  THL  POLLS  OH  AROUND  THE  VERNAL  EQUINOX 

2 IF ( UAbS  < ULC ) . LT • P 12  ) bO  TO  b 
RA=HA*Pi 

IF(UtC)  -5*3,4 

3 ULC=-ULC-PI 
bU  70  5 

4 ULt=Pl-UEC 

b IF (KA.GL.TPl ) K A=K A-TP 1 
iF(KA.LT.U.OO)  ka=ra+tpi 

return 

LNU 


15 


REFERENCES 


1.  K.  Haramundanis , Astron.  J.  72,  558  (1967). 

2.  F.  P.  Scott  and  C.  Smith,  "Conference  on  Photographic  Astrometric 
Technique,"  H.  Eichorn,  Editor,  NASA  CR-1825,  Washington,  D.C.  (1971), 
p.  167. 

3.  Smithsonian  Astrophysical  Observatory  (Staff  of).  Star  Catalog  (Publications 
of  the  Smithsonian  Institute  No.  4562,  Washington,  D.C.,  1966). 

4.  L.  G.  Taff,  "The  Handling  of  Star  Catalogs:  The  Transformations  of  Positions 

and  Proper  Motions,"  Technical  Note  1976-12,  Lincoln  Laboratory,  M.I.T. 

(17  February  1976). 

5.  L.  G.  Taff,  "Photometric  Reductions:  Theory  and  Practice,"  Technical  Note 

1976-35,  Lincoln  Laboratory,  M.I.T.  (17  December  1976). 


16 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (then  Doto  Entered ) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER 

2.  GOVT  ACCESSION  NO. 

3.  RECIPIENT'S  CATALOG  NUMBER 

ESD-TR-76-318 

4-  TITLE  (ond  Subtitle) 

5.  TYPE  OF  REPORT  4 PERIOD  COVERED 

The  Handling  and  Uses  of  the  SAO  Catalog 

Project  Report 

6.  PERFORMING  ORG.  REPORT  NUMBER 

Project  Report  ETS-6 

7.  AUTHOR^; 

«.  CONTRACT  OR  GRANT  numbers; 

Laurence  G.  Taff 

F19628-76-C-0002 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 
Lincoln  Laboratory,  M.I.T. 

P.  O.  Box  73 
Lexington,  MA  02173 

10.  PROGRAM  ELEMENT,  PROJECT.  TASK 
AREA  & WORK  UNIT  NUMBERS 

Program  Element  No.  63428  F 
Project  No.  2128 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

12.  REPORT  DATE 

Air  Force  Systems  Command,  USAF 
Andrews  AFB 
Washington,  DC  20331 

5 November  1976 

13.  NUMBER  OF  PAGES 

22 

14.  MONITORING  AGENCY  NAME  & ADDRESS  f'*/  different  from 

Controlling  Office) 

15.  SECURITY  CLASS,  (of  this  report ) 

Electronic  Systems  Division 
Han  scorn  AFB 

Unclassified 

Bedford,  MA  01731 

15a.  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 

16.  DISTRIBUTION  STATEMENT  (of  this  Report) 

Approved  for  public  release;  distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  (of  the  abstroct  entered  in  Block  20,  if  different  from  Report) 

18.  supplementary  notes 

None 

19.  KEY  WORDS  ( Continue  on  reverse  side  if  necessary  and  identify  by  block  number) 

artificial  satellites 
satellite  observation 

space  object  identification 

70.  ABSTRACT  ( Continue  on  reverse  side  if  necessary  ond  identify  by  block  number) 

The  complete  and  rigorous  update  procedures  applied  to  the  1950.0  epoch  SAO 
Catalog  are  described.  The  corrections  included  are  the  terms  of  elliptic  aberration, 
precession,  proper  motion,  and  foreshortening  (e.g.,  parallax  and  radial  velocity  terms). 
The  uses  of  the  updated  catalog  are  schematically  described  for  the  local  reduction  of  ar- 
tificial satellite  positions.  Finally,  an  annotated  FORTRAN  program  is  included. 

DD  1 1473  EDITION  OF  I NOV  65  IS  OBSOLETE 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (then  Dato  Entered) 


I 


\ 


